Restoring broken entanglement by injecting separable correlations 
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The distribution of entanglement is central in many protocols of quantum information and com- 
putation. However it is also known to be a very fragile process when loss and noise come into play. 
The inevitable interaction of the quantum systems with the external environment induces effects of 
decoherence which may be so strong to destroy any input entanglement, a phenomenon known as 
"entanglement breaking" . Here we study this catastrophic process in a correlated- noise environment 
showing how the presence of classical-type correlations can restore the distribution of entanglement. 
In particular, we consider a Gaussian environment whose thermal noise is strong enough to break 
the entanglement of two bosonic modes. In this scenario, we show that the injection of separable 
correlations from the same environment is able to reactivate the broken entanglement. This paradox- 
ical effect happens both in schemes of direct distribution, where a third party (Charlie) broadcasts 
entangled states to remote parties (Alice and Bob), and in schemes of indirect distribution which are 
based on the protocol of entanglement swapping, whose theory is here generalized. Furthermore, 
the amount of entanglement activated by the injection can be large enough to be distilled using 
one-way distillation protocols. As a result, entanglement distribution and its distillation are still 
possible in the presence of entanglement-breaking channels, as long as a sufficient amount of sepa- 
rable correlations can be identified in the environment. These findings pose fundamental questions 
about the intimate interplay between local and nonlocal correlations, and suggest new perspectives 
for quantum repeaters in the presence of memory channels and correlated-noise environments, such 
as those characterized by non-Markovian dynamics. 

PACS numbers: 03.65.Ud, 03.67.-a, 42.50.-p 



I. INTRODUCTION 

Entanglement is a fundamental physical resource in 
quantum information and computation [l|, Q ■ Once that 
two remote parties, say Alice and Bob, share a suitable 
amount of entanglement, they can implement a variety of 
powerful protocols, including the teleportation of quan- 
tum states [1, H[ and quantum gates [j| , as well as the 
distribution of unconditionally secure keys @, 0. The 
problem of entanglement distribution is therefore a cen- 
tral topic of investigation in the quantum information 
community. Unfortunately, this important resource is 
also fragile: Quantum systems inevitably interact with 
the external environment whose action may degrade or 
even destroy the presence of entanglement. Thus, in re- 
alistic implementations where the effect of decoherence is 
non-negligible, the distribution of entanglement becomes 
challenging and may need the help of distillation pro- 
tocols [g, l9( , where a large number of weakly-entangled 
states are transformed into a smaller number of strongly- 
entangled states via local operations and classical com- 
munications (LOCCs). 

Decoherence affects all the possible forms of entangle- 
ment distribution. In a scheme of direct distribution, 
there is a middle station (Charlie) possessing a bipartite 
system in an entangled state; one subsystem is sent to Al- 
ice and the other to Bob. Alternatively, we may consider 
a scheme of indirect distribution, known as entanglement 
swapping 

EMU, 

where the distribution is mediated by 
a measurement process. Here Alice and Bob each has 
a bipartite system prepared in an entangled state. One 
subsystem is retained while the other is sent to Charlie. 



At his station, Charlie detects the two incoming subsys- 
tems by performing a suitable measurement (known as 
the "Bell measurement" ) and communicates the classical 
outcome back to Alice and Bob. As a result of this pro- 
cess, the two subsystems retained by the remote parties 
are projected into an entangled state. 

In both the configurations, the distribution and dis- 
tillation of entanglement is possible as long as the deco- 
herent action of the environment is not strong. When 
decoherence is strong enough to destroy any input en- 
tanglement, the environment is described by the notion 
of entanglement breaking (EB) channel [III, Hi]]- A quan- 
tum channel £ is EB when its local action on one part 
of a bipartite state always results into a separable out- 
put state. In other words, given two systems, A and 
B, in an arbitrary bipartite state pab, the output state 
PAB' — (X4 ® £b)(pab) is always separable, where Ia 
is the identity channel applied to system A and £b is 
the EB channel applied to system B. Thus, if the input 
systems A and B were initially entangled (denoted by 
the notation A — B), the output systems A and B' are 
separable (denoted by the notation A\B'). 

Despite EB channels have been the subject of an in- 
tensive study by the community, they have been only 
analyzed under Markovian conditions of no memory for 
the environment. In other words, when the distribution 
involves two or more systems travelling in a completely 
decohering environment, these systems are typically as- 
sumed to be perturbed in an independent fashion, each 
of them subject to the same memoryless channel. 

For instance, consider the case of direct distribution 
depicted in Fig. HJl). In the standard memoryless de- 
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scription of the environment, the entangled state pab of 
the input systems A and B is subject to a tensor prod- 
uct of channels £a® £b- In this case, there is no way to 
distribute entanglement among any of the parties if both 
£a and £b are EB channels. Suppose that Charlie tries 
to share entanglement with one of the remote parties by 
sending one of the two systems while keeping the other 
(one-system transmission). For instance, Charlie may 
keep system A while transmitting system B to Bob. The 
action of Z4 <g) £ b destroys the initial entanglement, so 
that systems A (kept) and B' (transmitted) are separable 
(A\B'). Symmetrically, the action of £a ®Ib destroys 
the entanglement between system A' (transmitted) and 
system B (kept), i.e., we have A'\B. 

Now suppose that Charlie sends both his systems to 
Alice and Bob (two-system transmission). This strat- 
egy must necessarily fail under the assumption of no 
correlations in the environment. Since the joint action 
of the two EB channels is given by the tensor product 
£a® £b = {£a <& 1-b){Za ® £b) quantum entanglement 
must necessarily be destroyed. In other words, since we 
have one-system EB (A\B' and A'\B) then we must have 
two-system EB (A'\B'). 



Memoryless Correlated-Noise 
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FIG. 1: Direct and indirect schemes for entanglement dis- 
tribution in memoryless and correlated-noise environments. 
Charlie is the middle station transmitting to or receiving sys- 
tems from Alice (left station) and Bob (right station) . Ellipses 
represent entangled states, black-circles represent channels, 
and detectors are Bell measurements. (1) Direct distribution 
in a memoryless environment. If the individual channels are 
EB, then no distribution of entanglement is possile. We have 
A\B' and A'\B which implies A'\B' . (2) Entanglement swap- 
ping in a memoryless environment. We have that a\A' and 
B'\b implies a\b. (3) Direct distribution in a correlated en- 
vironment. Despite A\B' and A'\B we have that A' — B' is 
possible. Furthermore, this can be realized by a separable en- 
vironment. (4) Entanglement swapping in a correlated envi- 
ronment. Despite a\A' and B'\b we have that a — & is possible. 
Furthermore, this is realizable by a separable environment. 

The previous reasoning can be extended to the case of 



indirect distribution as shown in Fig.[lj2). Since the en- 
vironment is memoryless (£a®£b), we have that the ab- 
sence of entanglement before the Bell measurement (a\A' 
and B'\b) is a sufficient condition for the swapping pro- 
tocol to fail, i.e., the remote systems a and b remains 
separable (&|6). Similarly to the previous case, if one- 
system transmission does not distribute entanglement, 
then two-system transmission cannot lead to entangle- 
ment generation via the swapping protocol. 

In our paper, we show that the previous implications 
for direct and indirect distribution of entanglement are 
false in the presence of a correlated-noise environment, 
i.e., two-system transmission can successfully distribute 
entanglement despite one-system transmission is subject 
to EB. To the best of our knowledge, this is the first work 
to study the process of EB in the presence of correlated- 
noise environments, proving the advantages of the EB 
channels with memory. 

In our study, we will show that there are typical phys- 
ical conditions under which the environmental correla- 
tions are able to reactivate the distribution of entangle- 
ment, therefore "breaking entanglement-breaking" . The 
most remarkable finding is that we do not need to con- 
sider an entangled environment: The injection of separa- 
ble correlations from the environment is able to restore 
the distribution of entanglement. 

To better clarify these points, consider the schemes 
of direct and indirect distribution in the presence of 
a correlated-noise environment which are depicted in 
Figs. nj3) and (4). In the scheme direct distribution of 
Fig. [TJ3), an input entangled state pab is jointly trans- 
formed into an output state pa'B' = £ab(pab)- We as- 
sume that the dilation of the composite channel £ ab is 
realized by introducing a two-system environment, E\ 
and E2, in a bipartite state pe^e 2 -, which interacts with 
the systems via two identical unitaries Uae x (transform- 
ing A and E\) and Ube 2 (transforming B and E2). In 
other words, we have 

PA'B' = Tr El E 2 [(UAE 1 ®U B e 2 ) 

X {pAB ® PE X E 2 ){u AEl ®ul E2 )]. (1) 

If the environmental state is not tensor product, i.e., 
PEiE 2 7^ P-Ei ® Pe 2 , then the composite channel cannot 
be decomposed into memoryless channels, i.e., £ab 7^ 
£a ® £b- From the dilation given in Eq. ([T]), we can al- 
ways define the reduced channels, £a and £b, acting on 
the individual systems. For instance, if only system B is 
transmitted, then we have the evolved state 

PAB' = (1a ®£b)(pab) 

— ^e 2 [{Ia ® Ube 2 )pab ® Pe 2 {Ia®uI E2 )], (2) 

where pe 2 = Tte 1 (pe 1 e 2 )- A similar formula holds for 
the evolution of the first system p^B — (£a®^b)(pab)- 
In our work we show that, despite Ia®£b and £a®Ib 
are EB channels (so that A\B' and B\A'), the composite 
channel £ab is able to preserve entanglement (so that 
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A 1 — B' is possible). In other words, we show a para- 
doxical situation where Charlie is not able to share en- 
tanglement with Alice or Bob, but still he can distribute 
entanglement to them. This is clearly an effect of the in- 
jected correlations coming from the environmental state 
PEiE 2 - As mentioned earlier, our main finding is that 
these correlations do not need to be strong: Entangle- 
ment distribution can be activated by separable corre- 
lations, i.e., by an environment which is in a separable 
state PEiE 2 - 

This paradoxical result can also be extended to en- 
tanglement distillation, which typically requires stronger 
conditions than entanglement distribution (demanded by 
the existence of effective distillation protocols). Despite 
the individual channels are EB, their combination into a 
separable environment enables Charlie to distribute dis- 
tillable entanglement to Alice and Bob. This is easy to 
prove for an environment with finite memory, which can 
be decomposed as £ab ® Sab ® • • • 

In our investigation, we also consider the case of en- 
tanglement swapping in a correlated-noise environment 
as depicted in Fig. [IJ4). Here Alice and Bob have two 
entangled states, p a A and psb, respectively. Systems a 
and b are retained, while systems A and B are transmit- 
ted to Charlie, therefore undergoing the joint quantum 
channel Sab- Before the Bell measurement, the global 
state is described by 

PaA'B'b = (la ® £ AB ® T h )(p aA ® PBb) 

= Tr Bl£ ; 2 [U (paA ® PE X E 2 ® PBb)U^] , (3) 

where U = I a ® Uae x ® Ue 2 b ® h- As before, we con- 
sider the case where the reduced channels, £a and £b, 
are EB channels, so that no entanglement survives be- 
fore the Bell measurement (a\A' and B'\b). If the en- 
vironment has no memory (pEiE 2 = PEi <£> Pe 2 ) there 
is no way to distribute entanglement to Alice and Bob 
(a\b). By contrast, if the environment has memory 
{pExE 2 Pe x ® Pe 2 ), then entanglement distribution is 
possible (a — b) and this distribution can be activated by 
a separable environmental state PE\E 2 ■ Thus, we have 
the paradoxical situation where no bipartite entangle- 
ment survives at Charlie's station (a\A' and B'\b), but 
still the swapping protocol is able to generate remote en- 
tanglement at Alice's and Bob's stations (a—b) thanks to 
the separable correlations injected by the environment. 

As before, these separable correlations can be strong 
enough to distribute distillable entanglement to the re- 
mote parties. Note that this su gge sts a preferred mech- 
anism for quantum repeaters [15| in the presence of 
correlated-noise environments, where Charlie first per- 
forms the swapping protocol on each pair of systems re- 
ceived from the remote parties, followed by Alice and 
Bob performing entanglement distillation on the whole 
collection of swapped states. While this mechanism can 
be successful, the other standard approach of distillation 
followed by swapping cannot work since Charlie is not 
able to distill any entanglement with Alice or Bob. 



Our results are derived for continuous-variable sys- 
tems, i.e., quantum systems with an infinite-dimensional 
Hilbert space flml8l |. In particular, we have consid- 
ered the bosonic modes of the electromagnetic field (so 
that the various systems depicted in Fig. [T] can be in- 
terpreted as bosonic modes.) The input modes are pre- 
pared in Gaussian states with Einstein-Podolsky-Rosen 
(EPR) correlations (lg, |l9[ , which are the most typical 
form of continuous variable entanglement, and they are 
assumed to evolve under the action of a Gaussian envi- 
ronment. This type of environment is modelled by two 
beam splitters which mix the travelling modes, A and B, 
with two environmental modes, E\ and E2, prepared in a 
bipartite Gaussian state pe x e 2 (separable or entangled). 
The reduced channels, £a and £3, are two lossy chan- 
nels whose transmissivities and thermal noises are such 
to make them EB channels. The adoption of bosonic sys- 
tems and Gaussian channels is interesting not only for 
the potential practical applications, but also because the 
Gaussian formalism enables us to find closed analytical 
results. 

In terms of potential impact, our work opens new per- 
spectives for entanglement distribution and distillation 
in correlated-noise environments and memory channels, 
where the presence of correlations can be exploited to 
recover from the breaking of entanglement. Since the 
entanglement restoration can be achieved by the injec- 
tion of separable correlations, our work poses fundamen- 
tal questions on the intimate relations and joint dynam- 
ics of local and nonlocal correlations. It is important to 
note that memory channels and correlated (in particu- 
lar, non-Markovian) environments are present in a wide 
series of practical scenarios, including spin chains [20| . 
atoms in optical cavities [21], quantum dots in photonic 
crystals [22] , photonic propagation through linear optical 
systems [23H25| and atmospheric turbulence [2a - l28| |. 

The paper is structured as follows. In Sec. [II] we pro- 
vide basic notions on bosonic systems, Gaussian states 
and their main properties, which are useful for under- 
standing the remainder of the paper. (This section may 
be skipped by those readers who are familiar with contin- 
uous variable quantum information). In Sec. lIIII we intro- 
duce and characterize the basic model of correlated Gaus- 
sian environment, which directly generalizes the stan- 
dard model of thermal-loss environment. We identify the 
physical conditions under which the correlated Gaussian 
environment is separable or entangled, discussing spe- 
cific cases where the quadrature correlations are sym- 
metric or asymmetric. In Sec. IIV1 we study the direct 
distribution of entanglement in the presence of the corre- 
lated Gaussian environment and assuming the condition 
of one-system EB. We provide the regimes of parame- 
ters under which remote entanglement is activated by 
environmental correlations (in particular, separable cor- 
relations) and the stronger regimes where the generated 
remote entanglement is also distillable. The nature of the 
environmental correlations is also discussed in terms of 
quantum discord and the specific dynamics of the EPR 
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correlations is analyzed. In Sec.fV] we generalize the the- 
ory of entanglement swapping to the correlated Gaussian 
environment (with more details provided in appendix [A]) . 
We first consider the swapping of EPR correlations and 
then the swapping and distillation of entanglement, find- 
ing the regimes of parameters where swapping and distil- 
lation are successful despite the EB condition. Potential 
implications for quantum repeaters are also discussed. 
Finally, Sec. IVIl is for conclusion and discussion. 

II. BASIC NOTIONS ON BOSONIC SYSTEMS 
AND GAUSSIAN STATES 

The most important systems in continuous variable 
quantum information are the bosonic modes of the elec- 
tromagnetic field. A bosonic mode is a quantum system 
with an infinite-dimensional Hilbert space and described 
by a pair of quadrature operators: Position q and momen- 
tum p, satisfying the commutation relation [q,p] = 2i. 
More generally, a bosonic system of n modes is described 
by a vector of 2n quadrature operators 

x T := (qi,Pi, ■ ■ .,q n ,Pn) , ( 4 ) 

satisfying [xi, Xj] = 2iQij, where i,j = 1, . . . , 2n and f2y 
is the generic element of the symplectic form 

"■-e(-ii)- < 5 > 

fc=i v 7 

In experimental quantum optics, bosonic modes are 
typically prepared in Gaussian states. By definition, a 
quantum state p is "Gaussian" if its Wigner phase-space 
representation is Gaussian [l8j . This feature implies that 
any measurement of the quadratures provides outcomes 
which arc distributed according to a Gaussian statistics. 
From a theoretical point of view, a Gaussian state is very 
easy to characterize, being completely characterized by 
its first and second-order statistical moments. 

The first-order moment is known as the mean value, 
and is defined by x := (x), where (O) :— Tr(Op) denotes 
the average of the arbitrary operator O on the state p. 
The second-order moment is known as the covariance ma- 
trix (CM) V, with generic element 

V tJ :=±{{Axi,Axj}) , (6) 

where Aii := the deviation and { , } is the anti- 

commutator. Note the diagonal elements Vu correspond 
to the variances of the quadratures 

V( Xi ) := (Ax*) = (x*) - xl . (7) 

The CM is a 2n x 2n real symmetric matrix, which must 
satisfy the uncertainty principle [291 ] 

V + itl > , (8) 



implying the positive-definiteness V > [30]. 

The simplest examples of Gaussian states are ther- 
mal states. A single-mode thermal state has zero mean 
(x = 0) and CM equal to V = (2n + 1)1, where I is 
the 2x2 identity matrix and n is the average number 
of thermal photons. This state collapses to the vacuum 
state for n — 0. Multimode thermal states can be eas- 
ily constructed by tensor product. Remember that the 
operation of tensor product for the states p\ ® pi cor- 
responds to the operation of direct sum for the CMs 
Vi ffi V 2 . Then, the opposite operation of partial trace 
pi = Tr 2 (pi2) corresponds to collapsing the global CM 
V12 into the diagonal block associated with (qi,pi). 

In quantum optics' labs, Gaussian states are typically 
processed by Gaussian operations. The simplest ones are 
Gaussian unitaries, defined as those unitary operators 
which transform Gaussian states into Gaussian states. 
At the level of the statistical moments, the action of a 
Gaussian unitary p — > UpU* corresponds to 

x^ Sx + d, SVS T , (9) 

where d is a real displacement vector and S is a symplec- 
tic matrix, i.e., a real matrix preserving the symplectic 
form Sr2S T = f2. In the Heisenberg picture, a Gaussian 
unitary corresponds to the affine map 

x^Sx + d. (10) 

A Gaussian unitary is said to be a canonical unitary 
when d = 0. The most important example of canon- 
ical unitary is the beam-splitter transformation, which 
involves two bosonic modes. This unitary is character- 
ized by a transmissivity parameter < r < 1 and corre- 
sponds to the symplectic matrix 

Other important examples of Gaussian operations are 
the Gaussian channels [3l|. These channels describe all 
those cases where bosonic systems and environment in- 
teract by means of linear and/or bilinear Hamiltonians. 
The action of a Gaussian channel p — > £ (p) corresponds 
to the following transformation for the CM of the state 

V -> KVK T + N , (12) 

where K and N = N T are 2n x 2n real matrices, satis- 
fying the condition N + zf2 — iKfJK T >0||. 

The most important example of Gaussian channel is 
the lossy channel, which is typically used to model the op- 
tical propagation through dissipative linear media. When 
a single bosonic mode propagates through a lossy chan- 
nel, its CM is transformed by Eq. (TT2|) with 

K = y/?l, N= (l-r)(2n+l)I , (13) 

where < r < 1 is the transmissivity of the channel and 
n > its thermal number. A single-mode lossy channel 
can be dilated into a beam-splitter with transmissivity r, 
which mixes the incoming mode with an environmental 
thermal mode containing n mean photons. 
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A. Symplectic spectrum 

A central result in the theory of Gaussian states is 
Williamson's theorem [H, |33|. Given a CM V, there 
always exists a symplectic matrix S realizing the diago- 
nalization 



v = sws T , w = 0^ fe i , 



(14) 



fc=i 



where the diagonal matrix W is known as the Williamson 
normal form, and {vi,--- ,is n } are the n symplectic 
eigenvalues of the CM. Since det S = 1, we have that 



dctV =Y[v, 



(15) 



fe=i 



Using the symplectic spectrum, we can write the un- 
certainty principle of Eq. (|SJ in the equivalent form 

V > 0, v\ > 1 , (16) 

which clearly implies Vf- > 1 [34j . The symplectic spec- 
trum fully determines the entropic and purity properties 
of the Gaussian states. In fact, the von Neumann entropy 
S(p) := — Tr(plogp) of an arbitrary n-mode Gaussian 
state is expressed by [l8[ 



k=l 



where 



, . . i+l, x + 1 x — 1 x — 1 

Kx) ■= — =— lo § —5 tt- lo s —5— 



(17) 



(18) 



The basis of the logarithm can be taken to be 2 for bits or 
the Euler's number "e" for nats. The purity of a Gaussian 
state is given by 



1 



VdetV 



IK 1 



(19) 



fc=i 



so that the state is pure (mixed) iff dctV = 1 (> 1). 

In our paper, we will consider situations where the 
symplectic eigenvalues are large. In this case, it is useful 
to use the asymptotic expansion 



h(x) ~ log^x + O ( - 
2 \ x 



(20) 



which is valid for large x. For instance, if the whole 
symplectic spectrum is diverging [y^ large for any fc), 
then we can use the expansion of Eq. (|20|) and write the 
following asymptotic formula 

5(p)~log(-J VdetV = log — + nlog- , (21) 

where the entropy is simply related to the CM and the 
purity of the Gaussian state. 



B. Two- mode Gaussian states and bipartite 
entanglement 

Since we will consider the distribution of bipartite 
entanglement, we devote a specific section to Gaussian 
states of two bosonic modes and their separability prop- 
erties. Let us consider two modes, say A and B, with 
quadrature vector x T := (qA,PA,qB,PB)- We assume 
that these modes are described by a zero-mean Gaussian 
state pab- Since x = 0, this state is fully characterized 
by its CM, that we can always express in the blockform 



V = 



A C 
C T B 



(22) 



where A, B and C are 2x2 matrices. Finding the sym- 
plectic spectrum is straightforward, since [3, |36| 



'A± VA 2 -4dctV 



(23) 



where A := det A + dctB + 2dctC. The uncertainty 
principle is then equivalent to the bona-fide condition 



V > 0, 



> 1 



(24) 



where v- is the smallest symplectic eigenvalue. 

It is very easy to study the separability properties. It 
is sufficient to compute the smallest partially-transposed 
symplectic (PTS) eigenvalue e. This eigenvalue can be 
computed using the formula of u_, given in Eq. (|23p . 
where we replace A with A = det A + det B — 2 det C. 
Then, the Gaussian state is separable (entangled) if and 
only if e > 1 (e < 1). More strongly, this eigenvalue 
provides a quantification of entanglement, being an en- 
tanglement monotone for two-mode Gaussian states. In 
fact, it is monotonically related to the log-negativity [jjTj 



£ = max {0, — loge} 



(25) 



which is itself an entanglement monotone. 

The log-negativity £ also provides an upper-bound to 
the distillable entanglement, i.e., the average number 
of entanglement bits per copy which can be extracted 
from infinitely-many copies of the state pab <8> Pab ® • • ■ 
More interestingly, we can consider a lower-bound to the 
distillable entanglement. According to the hashing in- 
equality I38l . [39j , this is given by the coherent informa- 
tion da EJ 



I(A)B) = S(p B ) - S(pab) 



(26) 



where pb — Tta(pab) is the reduced state of Bob, here 
corresponding to a zero-mean Gaussian state with CM 
B. Using Eq. (fT7|) . the coherent information becomes 



I(A)B) = h{v B ) - h{vJ) - h(u+) 



(27) 



where vb = vdet B is the symplectic eigenvalue of B, 
and {j/_,z/ + } is the symplectic spectrum of V. In the 
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case where these eigenvalues are all very large, we can 
use the expansion of Eq. (I2U1) to get the formula 



I(A)B)~\og- 



detB 
detV 



(28) 



Remember that the coherent information is a lower 
bound to the distillable entanglement that can be ob- 
tained by means of one-way protocols between Alice and 
Bob. In these protocols, Alice applies a quantum in- 
strument to her copies, i.e., a collection of completely- 
positive trace-preserving maps labelled by a classical in- 
dex k. Then, she classically communicates k to Bob, who 
applies a conditional quantum operation on his copies 
(see Ref. [38[ for more details). 



C. EPR correlations 

The most important example of two-mode Gaussian 
state is the two-mode squeezed vacuum state, also known 
as the EPR state. This state represents the most common 
source of continuous-variable entanglement, and will be 
adopted as such a source in our study of entanglement 
distribution. 

An EPR state has zero mean and CM given by 



V = 



pi p'Z 

p'z pi 



where p > 1, pi := sj p 2 — 1, and 



(29) 



(30) 



The quadratures x T := (qA,PA,<lB,PB) have equal vari- 
ances p and maximal correlations (qAQB) = — {paPb) = 
p' . The parameter p can be used to quantify the entan- 
glement between the two modes A and B. In fact, the 
smallest PTS eigenvalue takes the form 



V2/x (a* - /x') - 1 



(31) 



which is monotonically decreasing in p. We can easily 
check that we have entanglement (e < 1) for any p> 1. 

The EPR correlations can be characterized in terms of 
the variance of the EPR quadratures 



QA ~ Qb „ Pa + Pb 
c i- '■= l=— ' P+ z 



V2 



V2 



(32) 



The EPR condition corresponds to V(q-) — V(p+) < 1, 
which means that the quadrature correlations between 
the two modes are below the vacuum noise. For the EPR 
state described by Eq. (|2"9"|) indeed we have 



(33) 



which is less than 1 for any p > 1. For this type of state, 
EPR correlations and entanglement are equivalent condi- 
tions. In general, the presence of EPR correlations in an 



arbitrary two-mode Gaussian state represents a sufficient 
condition for entanglement. 

In the limit of large entanglement (p — > oo), the EPR 
state becomes ideal. In fact, we have V(q~) — V(p+) — > 
0, which corresponds to realizing the ideal EPR corre- 
lations qA = q_B and pa = —Pb, where positions (mo- 
menta) are perfectly correlated (anticorrelated). Note 
that an alternative EPR state has — Z in the place of Z 
in Eq. (|2"9"|) . This other state has EPR correlations in the 
quadratures <?+ and p- , therefore tending to realize the 
ideal conditions qA — —<1b and pa = Pb- In our work, we 
do not consider this alternative state, but similar results 
can easily be derived. 



D. Gaussian quantum discord 

Quantum entanglement is synonymous of quantum 
correlations for pure states. For mixed states the scenario 
is more subtle. In fact, separable mixed states may still 
contain features of quantumness, one of which is known 
as quantum discord [H, |43j]. Quantum discord derives 
from the disagreement between the quantum generaliza- 
tion of two information quantities. 

Classically, we know that the mutual information of 
two random variables A and B can be defined as 

I(A : B) = H(A) + H(B) — H(A, B) , (34) 

where H is the Shannon entropy. Equivalently, using the 
chain rules for the entropy, we can write 



I {A : B) = H(A) - H(A\B) 



(35) 



where H(..\..) is the conditional Shannon entropy. 

In quantum mechanics, the previous quantities do not 
have equivalent generalizations. Given a bipartite quan- 
tum state pab, its quantum mutual information is ex- 
pressed in terms of the von Neumann entropy as 



AB 



S(p A ) + S(p B ) - S(pab) 



(36) 



This is the generalization of Eq. (f3"4"]) and accounts for all 
the correlations present in the quantum state. Different 
is the generalization of Eq. (|35|) which is given by 

C{ PAB ) = S(p A ) - MH M {A\B) , (37) 

where M = {Mk} is a POVM acting on system B and 

H M (A\B) :=5^j» fc S(pA|fc) , (38) 
k 

where pk is the probability of the outcome k and pA\k is 
the conditional state of A. The quantity C(pab) quan- 
tifies the classical correlations in the state correspond- 
ing to the amount of common randomness which can 
be extracted using one-way classical communication pij . 



Quantum discord is defined as the difference between the 
total correlations and the classical correlations, i.e. 

D{ PAB ) = I( PAB ) - C{ PAB ) ■ (39) 

Note that Eq. (|3T[) is generally non-symmetric under 
system permutation. As a result, we have that the B- 
type quantities C(pba) and D(pba) are generally differ- 
ent from the A- type quantities of Eqs. (|3"T)) and (|3"9"]l . 

For bosonic Gaussian states, we can define the so- 
called Gaussian quantum discord where the minimization 
is restricted to Gaussian POVMs. Consider a two-mode 
Gaussian state with CM in the blockform of Eq. (f22j). 
Then, classical correlations and Gaussian quantum dis- 
cord are given by the formulas djj |4(| 

C(pab) = h (Vdct A) - E , (40) 

and 

D(pab) = h (VdetB) - h(vJ) - h(u+) + E , (41) 

where the expression of the term E is given in Ref . [46j . In 
our paper, we consider symmetric Gaussian states (A = 
B) so that there is no ambiguity in which type of discord 
and classical correlations we are using. 

III. GAUSSIAN ENVIRONMENT WITH 
CORRELATED NOISE 

In this section, we introduce and fully characterize our 
model of correlated Gaussian environment, which repre- 
sents the simplest and most direct generalization of the 
standard memoryless Gaussian environment with losses 
and thermal noise. 

Let us consider a Gaussian decoherence process which 
affects two bosonic modes, A and B, in a symmetric way. 
In order to describe the introduction of losses and ther- 
mal noise, we consider two beam splitters (with the same 
transmissivity r) which combine modes A and B with two 
environmental modes, E\ and i?2, respectively. These an- 
cillary modes are prepared in a zero-mean Gaussian state 
PE t E 2 which is symmetric under E\-Ei permutation. In 
the standard memoryless model depicted in Fig. (3Ji) , the 
environmental state is taken to be in a tensor-product 
PEiE 2 — P ® Pi meaning that E\ and E2 are fully in- 
dependent. In particular, p is a thermal state with CM 
col, where the noise variance lo — 2n + 1 quantifies the 
mean number of thermal photons h entering the beam 
splitter. Each interaction is then equivalent to a lossy 
channel with transmissivity r and thermal noise u). 

In our work, we generalize this Gaussian process to 
include the presence of correlations between the environ- 
mental modes as depicted in Fig. HJii). The simplest 
extension of the model consists of taking the ancillary 
modes, E\ and E2, in a zero-mean Gaussian state pe^e 2 
with CM given by the symmetric normal form 

V ElE2 (u,g,g')= ( £ ^) , (42) 




FIG. 2: Models of Gaussian environments, (i) Memoryless 
Gaussian environment with losses r and thermal noise w. 
(ii) Correlated Gaussian environment, with losses r, thermal 
noise lu and correlations G. 



where lo > 1 is the thermal noise variance associated with 
each ancilla, and the off-diagonal block 

G = (%,), (43) 

accounts for their correlations. 

Clearly, when we consider the two interactions A — E\ 
and B — E2 separately, the environmental correlations 
are washed away. In fact, by tracing out E2, we are left 
with mode E\ in a thermal state (Ve x = loT) which is 
combined with mode A via the beam-splitter. In other 
words, we have again a lossy channel with transmissivity 
r and thermal noise u). The scenario is identical for the 
other mode B when we trace out E%. However, when we 
consider the joint action of the two environmental modes, 
the correlation block G comes into play and the global 
dynamics of the two travelling modes becomes completely 
different from the standard memoryless scenario. 

Before studying the system dynamics and the corre- 
sponding evolution of entanglement, we need to charac- 
terize the correlation block G more precisely. In fact, 
the two correlation parameters, g and <?', cannot be 
completely arbitrary but must satisfy specific physical 
constraints. These parameters must vary within ranges 
which make the CM of Eq. (l4"2l a bona- fide quantum CM. 
Given an arbitrary value of the thermal noise U) > 1, the 
correlation parameters must satisfy the following three 
bona-fide conditions 

\g\<U, \g'\<U), u 2 +gg'-l>u;\g + g'\. (44) 

The proof of Eq. (|4"H) is easy. We impose Eq. to 
Ve 1 e 2 ( uj , g, g'). The positivity Ve ± e 2 > is equivalent 
to the positivity of the principal minors of the matrix. 
The positivity of the first two minors is trivially implied 
by cj > 1. The third minor gives 

(u g\ 

det OwO >0» u{u? - g 2 ) > 0, (45) 

\g lo ) 

which is equivalent to 

\g\ < <* ■ (46) 
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The fourth minor corresponds to the determinant 

dctV ElE2 = {uj 2 - g 2 )(uj 2 - g' 2 ) , (47) 

and its positivity det Vg^a > leads to the condition 

\g'\ < uj . (48) 

Finally, using A = 2(uj 2 +gg') and Eq. (g7J|, we find that 
v 2 _ > 1 is equivalent to [47j | 

w 2 +.g.g'-l >w|.g + .g'| . (49) 

In conclusion, an environmental CM of the form (|42"j) 
with thermal noise uj > 1 is a bona-fide CM if the cor- 
relation parameters g and satisfy the three bona-fide 
conditions of Eq. flUJ. Using Eqs. §T§§ and (|IT|). we also 
see that its purity is given by 

mB2 = [( w 2 - 5 2 )( w 2 - 5 ' 2 )]- 1/2 . (50) 

A. Separability properties 

Once that we have fully clarified the bona-fide condi- 
tions for the environment, the next step is to characterize 
its separability properties. 

For this aim, we compute the smallest PTS eigenvalue 
associated with Ve 1 e 2 - After simple algebra, we get 

£= V^ 2 - gg' - uj\g - g'\ . (51) 

Provided that Eq. ([4"4"]l is satisfied, the separability con- 
dition e > 1 is equivalent to 

uj 2 -gg'-l >io\g-g'\ . (52) 

The various conditions of bona-fide and separabil- 
ity can be combined together. An environment of the 
form (I42[) with thermal noise uj > 1 is bona-fide and sep- 
arable when the correlation parameters satisfy 

\g\ < uj, \g'\ <uj, to 2 - 1 > max{r_, r + } . (53) 

where 

r_ ~u\g + g'\-gg' , (54) 

r+ :=uj\g-g'\+gg' . (55) 

By contrast, it is bona-fide and entangled when 

\g\ < u, \g'\ <u, r_ < uj 2 - 1< T+ . (56) 

To better clarify the structure of the environment, we 
provide a numerical example in Fig. [3J In this figure, 
we consider the correlation plan which is spanned by the 
two parameters g and g' . For a given value of the thermal 
noise uj, we identify the subset of points which satisfy the 
bona-fide conditions of Eq. flU} . This subset corresponds 
to the white area in the figure. Within this area, we then 
characterize the regions which correspond to separable 
environment (area labelled by S) and entangled environ- 
ment (areas labelled by E). 




-2-1 1 2-2-1 1 2 



9 9 

FIG. 3: Left. Correlation plan (g, g') for the Gaussian envi- 
ronment, corresponding to thermal noise uj — 2. The black 
area identifies forbidden environments (correlations are too 
strong to be compatible with quantum mechanics). White 
area identifies physical environments, i.e., the subset of points 
which satisfy the bona-fide conditions of Eq. (|44[) . Whitin 
this area, the inner region labbeled by S identifies separable 
environments [Eq. (|53[) is satisfied] while the two outer re- 
gions identify entangled environments [Eq. (|56p is satisfied]. 
Right. Correlation plan displaying the SC environments (bi- 
sector g' = g) and the AC environments (bisector g = —g). 
The MSC and EPR environments are also diplayed (extremal 
points on the bisectors). 



B. Specific types of environments 

Here we specify our model of Gaussian environment 
to particular cases with specific correlation properties. 
These types of environment will be taken as examples 
for deriving analytical results. 

Before discussing these cases, it is trivial to say that, 
for g' = g = 0, our environment collapses into a memo- 
ryless Gaussian environment, described by two identical 
lossy channels with transmissivity r and thermal noise 
uj. This is represented by the origin in Fig. [3] Also 
note that the condition uj = 1 corresponding to vacuum 
noise is only compatible with this environment. In other 
words, for uj = 1, we must necessarily have g' = g = 0, 
since this is the unique solution which is compatible with 
the bona-fide conditions of Eq. (|44| . 

The first type of correlated environment, that we call 
"symmetrically correlated" (SC) corresponds to taking 
g' = g. This means that positions and momenta of the 
two environmental modes are correlated in exactly the 
same way. On the correlation plan, this environment is 
represented by one of the points lying on the bisector of 
the first and third quadrants (see Fig. [3]). 

It is easy to check that the bona-fide conditions of 
Eq. (|44[) simplify to the unique inequality 

\g\ < w - 1 , (57) 

so that the maximal symmetrical correlations (MSC) are 
given by g = uj — 1 oi j = 1 - uj, which are the two 
extremal points indicated in Fig. |3J As we can see from 
the figure, the SC environments are always separable, so 
that Ei and are correlated but not entangled. In 
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fact, the separability condition of Eq 
yw 



521) becomes \g\ < 
2 — 1 which is always satisfied by uj > 1 and Eq. 



The SC environment is separable but generally mixed, 
since Eq. (|50|) becomes pe 1 e 2 < (2w — which is less 
than 1 for any oj > 1 . It is pure only for u> = 1 when it col- 
lapses into a memory less environment. Using Eqs. (|4"0")) 
and pTj) . we investigate the nature of its separable cor- 
relations in terms of classical correlations and quantum 
discord. The amount of classical correlations is equal to 



C = h(u) 



9 



1 



while the Gaussian quantum discord is given by 

n 2 

D = h{uj) — h(u) — g) — h(u + g) + h 
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LO+1 



(58) 



(59) 



The second type of environment, that we call 
"asymmetrically correlated" (AC) corresponds to the 
condition g' = —g. This means that positions and mo- 
menta have opposite correlations, i.e., if positions are cor- 
related (anticorrelated), then momenta are anticorrelated 
(correlated). On the correlation plan, this environment 
is represented by one of the points lying on the bisector 
of the second and fourth quadrants (see Fig. [3]). 

In this case, the bona-fide conditions of Eq. (|4"4"|) sim- 
plify to the inequality 



\g\ < 



i 



(60) 



For maximal correlations |g| = \Jlo 2 — 1 we have an EPR 
state, which is pure and maximally entangled. Depend- 
ing on the sign of <?, we have two different EPR environ- 
ments: The positive EPR environment (g = y/co 2 — 1) 
with positions correlated, and the negative EPR envi- 
ronment (g = —\/u 2 — 1) with positions anticorrelated. 

In general, the AC environment can be separable or 
entangled, as evident from Fig. [3] By using Eq. (I5"2"j) . we 
see that it is separable for 

\g\<u-i, (61) 

while it is entangled for stronger correlations 

Lu-K\g\< \]io 2 - 1 . (62) 

Finally, we can compute the amount of classical correla- 
tions and Gaussian quantum discord. These quantities 
are respectively given by Eq. (|58[) and 



D = h{uj) - 2h ( y/u 2 - g 2 + h lu) - 
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LO+1 



(63) 



IV. 



DIRECT DISTRIBUTION OF 
ENTANGLEMENT 



Let us study the system dynamics and the entangle- 
ment propagation in the presence of a correlated Gaus- 
sian environment. Suppose that Charlie has an entan- 
glement source described by an EPR state pab with CM 



V(/i) given in Eq. (l29l) . We may consider the different 
scenarios depicted in Fig. @] Charlie may attempt to 
distribute entanglement to Alice and Bob as shown in 
Fig. E0), or he may try to share entanglement with one 
of the remote parties, as shown in Figs.[4|ii) and (hi). 



(i) 


Charlie 












o\ /u>\ 


A / 






E, E 2 \ 


Alice 


Bob 




FIG. 4: Scenarios for direct distribution of entanglement, (i) 
Charlie has two modes A and B prepared in an EPR state 
Pab- In order to distribute entanglement to the remote par- 
ties, Charlie transmits the two modes through the correlated 
Gaussian environment characterized by transmissivity r, ther- 
mal noise uj and correlations G. (ii) Charlie aims to share en- 
tanglement with Alice. He then keeps mode B while sending 
mode A to Alice through the lossy channel £a- (iii) Charlie 
aims to share entanglement with Bob. He then keeps mode A 
while sending mode B to Bob through the lossy channel £b- 

Let us start considering the scenario where Charlie 
aims to share entanglement with one of the remote par- 
ties (one- mode transmission). In particular, suppose that 
Charlie wants to share entanglement with Bob (by sym- 
metry the derivation is the same if we consider Alice). 
For sharing entanglement, Charlie keeps mode A while 
sending mode B to Bob as shown in Fig. @J hi). The ac- 
tion of the environment is therefore reduced to I a ® £b, 
where £b is a lossy channel applied to mode B. It is easy 
to check that the output state pab' shared by Charlie and 
Bob is Gaussian with zero mean and CM 



AB' 



where 



/ii p'y/TZ \ 

fi'y/TZ Xl J 



Tfl + (1 - T)UJ 



(64) 



(65) 



In fact, the global input state pab <£> Pe 2 nas CM 
V(/x)©ojI. This CM is subject to the symplectic transfor- 
mation I © S(r), where S(r) is the beam splitter matrix 
of Eq. (fTTj) acting on modes B and E% . Tracing out the 
environment from the output CM, we get Eq. (|64l) . 

Remarkably, we can compute closed analytical formu- 
las in the limit of large p, i.e., large input entanglement. 
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In this case, the entanglement of the output state is quan- 
tified by the PTS eigenvalue 



1 



1 



(66) 



The EB condition corresponds to e < 1, which provides 

(67) 



1+T 

uj > - := w E B , 



or equivalcntly 



1 -r 



n > 



1 



(68) 



Despite the EB condition of Eq. (1571) has been derived 
for an EPR input, it is actually completely general. In 
other words, a lossy channel £b with transmissivity r 
and thermal noise w > web destroys the entanglement 
whatever Charlie's input state pab is. In fact, Eq. (|6"Tj) 
corresponds exactly to the general EB condition for lossy 
channels which has been derived in Ref. 14|. It is there- 
fore clear that the threshold condition uj = web guaran- 
tees one-mode EB, i.e., the impossibility for Charlie to 
share entanglement with Alice or Bob. 

Now the central question is the following: Suppose that 
Charlie cannot share entanglement with the remote par- 
ties (one- mode EB), can Charlie still distribute entan- 
glement to them? In other words, suppose that the cor- 
related Gaussian environment has transmissivity r and 
thermal noise uj = web, so that the lossy channels Ea 
and £b are EB. Is it still possible to use the joint chan- 
nel Sab to distribute entanglement to Alice and Bob? In 
the following, we explicitly answer to this question. Fur- 
thermore, we will show that the distributed entanglement 
can also be distilled by one-way protocols. 

Let us derive the general evolution of the two modes A 
and B under the action of the environment [see Fig. \MX)} ■ 
Since the input EPR state pab is Gaussian and the envi- 
ronment is Gaussian, the output state pa> B' is also Gaus- 
sian. This state has zero mean and CM given by 



Va'B' = tVab + (1 - t)V EiE2 = 



xl H 
H xl 



(69) 



where 

H := r//Z + (1 — t)G . (70) 
For large p, we can easily derive the symplcctic spectrum 
u ± = ^(2uj + g> -g±\g + gi\)(l-T)Tp , (71) 



and the smallest PTS eigenvalue 



(l-T)y/(w-g)(LJ + g') 



(72) 



which quantifies the amount of entanglement distributed 
to Alice and Bob. In the same limit, we can easily com- 
pute the coherent information between the two remote 



parties. In fact, using the formula of Eq. (f28|) for diverg- 
ing spectra, we get 



I(A)B) = log 



1 



ee 



(73) 



Thus, remote entanglement is distributed for e < 1 and 
this entanglement is distillable for e < e _1 . 

Suppose that the environment has thermal noise to = 
Web (one-mode EB). Then, we can write 



e = y/[\ + T-{l- r).g][l + r + (1 - r)g'] 
~<r,g,g') (74) 

Answering our previous question corresponds to checking 
the existence of environmental parameters r, g and g' , for 
which e < 1. For a given value of the transmissivity r, 
we look for regions in the correlation plan {g,g') where 
remote entanglement is distributed (e < 1) and possibly 
distillable (e < e _1 ). This is done in Fig. [S] for several 
numerical values of the transmissivity. 

In Fig. [5J the environments identified by the gray re- 
gion, that we call the "activation area", allow Charlie 
to distribute entanglement to Alice and Bob (e < 1), 
despite it is impossible for him to share entanglement 
with the remote parties. In other words, these environ- 
ments are two-mode entanglement preserving (EP), de- 
spite they are one-mode EB. Furthermore, we can also 
identify environments which are able to generate distill- 
able entanglement (e < e _1 ). 

The most remarkable feature in Fig.[5]is represented by 
the presence of separable environments in the activation 
area. In other words, there are separable environments 
which contain enough correlations to activate the distri- 
bution of entanglement to Alice and Bob. Furthermore, 
for sufficiently high transmissivities and correlations, this 
entanglement can also be distilled by means of one-way 
protocols. Also note from Fig. [5] that the weight of sep- 
arable environments in the activation area increases for 
increasing transmissivities, with the entangled environ- 
ments almost disappearing for r = 0.9. 

It is important to stress that achieving the simultane- 
ous conditions of one-mode EB and two-mode EP is not 
surprising if we consider entangled environments. For 
instance, we may consider two beam-splitters with zero 
transmissivity, so that Charlie's state is completely re- 
flected into the environment and the (entangled) state 
of the environment is reflected to Alice and Bob. This 
scenario is certainly one-mode EB since Charlie has no 
chance of sharing part of his state with Alice or Bob, and 
it is also two-mode EP since the loss of Charlie's initial 
entanglement is replaced by the injection of entanglement 
from the environment, which is then distributed to Alice 
and Bob. 

By contrast, achieving the conditions of one-mode EB 
and two-mode EP with separable environments is sur- 
prising, because no injection of entanglement is present 
but still the environmental correlations are strong enough 
to restore the broken entanglement. 
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FIG. 5: Analysis of the remote entanglement on the correla- 
tion plan for different values of the transmissivity r = 0.3, 0.5, 
0.75, and 0.9 (from top left to bottom right). Corresponding 
values of the thermal noise are determined by the one-mode 
EB condition lj = ojeb- In each inset, the non-black area 
identifies the set of physical environments, which are divided 
into separable (S) and entangled (E) environments. The gray 
region is the activation area and identifies those environments 
which enable Charlie to distribute entanglement to Alice and 
Bob (e < 1). Within the activation area, the environments 
below the dashed curve enable Charlie to distribute distillable 
entanglement to the remote parties (e < e _1 ). 



A. Direct distribution in specifically correlated 
environments 

Here we analyze the scheme of direct distribution in 
the presence of the specific environments discussed in 
Sec. lIII Bl i.e., the SC and AC environments, correspond- 
ing to the bisectors of the correlation plan. We consider 
these environments at the EB threshold u> — web , so that 
they are one-mode EB. In the limit of large /i, we derive 
the regimes of parameters for which remote entanglement 
can be distributed and distilled. 

Let us start with the SC environment (g' = g) which is 
always separable. In this case, the remote entanglement 
of Alice and Bob is quantified by e — e(t, g, g) according 
to Eq. ([73]). This quantity takes the optimal value 



£msc - V(1-t)(1+3t) 



(75) 



when the environment is MSC, i.e., it has maximal cor- 
relation parameter 



|<?msc| = web — 1 



2r 

1 -T 



(76) 



From Eq. (|75j) we can easily check that remote entangle- 
ment is generated (e < 1) only if r > 2/3, and remote 
entanglement can be distilled (e < e^ 1 ) only if r > 0.96. 

For an arbitrary SC environment, remote entangle- 
ment is generated for r > 2/3 and 

IstegI < \g\ < |5msc| , (77) 
where the entanglement generation bound in given by 



\9EG\ 



Vt(t + 2) 

1 - T 



(78) 



Entanglement is then distillable when r > 0.96 and 

|pbd| < \g\ < IsmscI ■ 



(79) 



where the entanglement distillation bound is equal to 



|SEd| 



Ve 2 (l + r) 2 -l 



e(l-r) 

See Fig. [5] for a pictorial representation 



> |0eg| 



(80) 
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FIG. 6: Schematic representation showing the separability 
properties of the remote state in terms of the correlation pa- 
rameter of the SC environment. The environment is separable 
and satisfies the one-mode EB condition uj — ujeb- Here we 
consider high transmissivity r > 0.96 so that remote entan- 
glement can also be distilled. 

Now consider the AC environment (g' — —g) which is 
separable for \g\ < \g$\ and entangled for \g$\ < \g\ < 
|Sepr|, where 



\98\ 



web — 1, |5epr| 



EB 



Vt 

1 -T 



Alice and Bob's remote entanglement is quantified by 



£ = e(r, g, -g) = l + r-(l- r)g 



(82) 



which can be less than one only for positive values of the 
correlation parameter g. This asymmetry, which is also 
evident from Fig. [5l depends on the fact that Charlie's 
input state has EPR correlations of the type qa = <1b 
and pa — —f>B- These correlations tend to be preserved 
by AC environments with positive g (having correlations 
of the same type) while they tend to be destroyed by AC 
environments with negative g (having correlations of the 
opposite type). 

The optimal distribution of entanglement is clearly 
achieved for the positive EPR environment g — |<7epr|, 
for which we have £epr = (1 — V^") 2 - The optimal dis- 
tribution in a separable AC environment is achieved at 
the border value g = \g§\ where es = 1 — t. From the 
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expression of Eq. (|82|) , we see that remote entanglement 
is generated at any transmissivity r for values of the cor- 
relation parameter 

ffEG < 9 < |5epr| , (83) 
where the entanglement generation bound is equal to 
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(84) 



i 1 1 1 — ► g 

9eg 1 9s I |9eprI 

FIG. 7: Schematic representation showing the separability 
properties of the remote state versus those of the AC environ- 
ment in terms of the correlation parameter. The environment 
satisfies the one-mode EB condition oj — cjeb- 
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FIG. 9: Remote entanglement e, Gaussian quantum discord 
D and purely-classical correlations C are plotted as functions 
of the correlation parameter g, for the SC environment (left) 
and the AC environment (right). Transmission is r = 0.9 
and thermal noise is ujeb = 19 (one-mode EB). The SC en- 
vironment is always separable while the AC environment is 
separable for \g\ < 18 and entangled for 18 < \g\ < 18.97. 
Remote entanglement is generated (e < 1) when C is appre- 
ciably different from zero. Remote entanglement is optimal 
for the MSC environments in the left panel, and the positive 
EPR environment in the right panel. 



As shown in Fig. [7] the generation of remote entan- 
glement is possible in the presence of separable environ- 
ments (<?eg < I <7s I)- Entanglement distillation is possible 
for r > e _1 (\/e — l) 2 ~ 0.15 and correlation values 

<?ed < 9 < |5epr| , (85) 

where the entanglement distillation bound is given by 

1 + r- - 1 



.9ed 



1 — T 



> .9eg 



(86) 



At higher transmissivities r > 1 — e _1 ~ 0.63 we have 
that distillation is possible even for separable environ- 
ments (#ed < IssD- For a schematic see Fig. [5] 
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FIG. 8: Schematic representation showing the separability 
properties of the remote state versus those of the AC envi- 
ronment in terms of the correlation parameter. The environ- 
ment satisfies the one- mode EB condition oj — ojeb- Here we 
consider high transmissivity r > 0.63 so that remote entan- 
glement is distillable even in the presence of separable envi- 
ronments. 

Finally, it is interesting to investigate what kind of 
correlations are present in separable environments which 
are one-mode EB and two-mode EP. This is done in Fig. [5] 
for the SC and AC environments. Remote entanglement 
is generated (e < 1) when purely-classical correlations C 
are appreciably different from zero. The role of Gaussian 
quantum discord D seems to become more important for 
high values of g, in particular for the AC environment, 
where D approaches C when the environment becomes 
entangled, reaching the extreme value D = C for the 
EPR environments. 



B. Evolution of the EPR correlations 

In order to give another point of view to the dynamics 
of the process, we describe here the evolution of the EPR 
correlations under the effect of the correlated Gaussian 
environment [see Fig.[5Ji)]. We can easily write the input- 
output Bogoliubov transformations 



(87) 



X.4' = VrXA + n/T — TXEi 



where x = (q,p) T is a vector of quadratures. From these 
equations, we can extract the output EPR operators 



<L4' - IB' PA' + PB' 

1- '■= ^ . P+ = 



V2 



v2 



(89) 



Now, using the CM of the input EPR state and that of 
the environment, we can compute the variances 

' I V(p' + ) ) 

= r(/i-p')I + (l-r)M-ZG). (90) 
In the limit of [i 3> 1 , we have 

A — s- Aoo = (1 — t){u1 — ZG) , (91) 
and assuming the EB condition lj — web we get 

A 00 , EB = (l + r)I-(l-r)ZG . (92) 



Analyzing Eq. (|92j) we may guess the possibility of real- 
izing the EPR condition Aq^eb < I by adopting suitable 
choices of the correlation block G. 

Let us explicitly compare the different types of envi- 
ronments. For the memoryless environment (G = 0) we 
have Aqq^eb = (1 + t)I which means that the EPR vari- 
ances are always greater than or equal to one, i.e., EPR 
correlations do not survive. 



For the AC environment (G = gZ) we have Aqo.eb = 
[(1 + t) - (1 - r)g]I. It is easy to check that the EPR 
condition Aoo.eb < I is achieved by physical values of 
9 > 3eGj where <7eg is given in Eq. (I84[) . This means that 
the remote entanglement generated by this environment 
is always in the form of EPR correlations (of the same 
type of the original EPR state at Charlie's station). 

For the SC environment (G = gl) we have Aq^eb — 
(1 + t)I — (1 — r)gZ and we can check that the condi- 
tion Aoo^eb < I is not realizable by any choice of g. In 
particular, for the MSC environment g = cjeb — 1 we 
have 

Aoo,eb= ( 1_T 1 + 3r ) • (93) 

Thus, we see that the initial EPR correlations do not 
survive in this case. Nevertheless remote entanglement 
can be distributed in the presence of this environment. 



V. INDIRECT DISTRIBUTION OF 
ENTANGLEMENT 

In this section we consider the indirect distribution of 
entanglement, i.e., the protocol of entanglement swap- 
ping. We start with a brief review of this protocol in 
the ideal case of no noise. Then, we generalize its theory 
to the case of correlated-noise Gaussian environments, 
where we prove how entanglement swapping can be acti- 
vated in the presence of one-mode EB. 



A. Entanglement swapping in the absence of noise 

Consider two remote parties, Alice and Bob, who pos- 
sess two identical EPR states with CM given in Eq. (f2!)f . 
At Alice's station, the EPR state describes modes a and 
A, while at Bob's station it describes modes b and B. 
Alice and Bob keep modes a and 6, while sending modes 
A and B to Charlie, where a Bell measurement is per- 
formed. This means that the travelling modes A and B 
are combined in a balanced beam splitter whose output 
modes "— " and "+" are homodyned, with mode "— " 
measured in the position quadrature and mode "+" in 
the momentum quadrature. In other words, Charlie mea- 
sures the two EPR quadratures q- and p+ which are de- 
fined by Eq. (13"21) . The Bell measurement provides two 
outcomes, g_ and p+, which can be compacted into a 
single complex variable 7 := q_ + ip+. The classical vari- 
able 7 is finally communicated to Alice and Bob, with 
the result of projecting their remote modes a and b into 
a conditional state p a b\-y (see Fig. ITU|). 

Since the input states are pure Gaussian and the Bell 
measurement is a Gaussian measurement which projects 
pure states into pure states, we have that the remote con- 
ditional state p a b\j turns out to be a pure Gaussian state. 
This state has a measurement-dependent mean x = x(7) 



13 




FIG. 10: Entanglement swapping in the absence of noise. See 
text for explanation. 

which can always be deleted by Alice and Bob via condi- 
tional displacements. It is clear that these local unitaries 
do not alter the amount of entanglement in the state, 
as long as they are perfectly implemented (see Ref. [48| 
for a general analysis which includes imperfect displace- 
ments). The conditional CM "V a b\-y can be computed us- 
ing the input-output formula for Gaussian entanglement 
swapping which has been proven in Ref. [H[ . We get 




Its smallest PTS eigenvalue is equal to e — jj," , which 
means that remote entanglement is always generated for 
entangled inputs (p > 1). Furthermore, remote en- 
tanglement is present in the form of EPR correlations 
since the two EPR quadratures q- = {q a — qb)/V% and 
p+ = {p a + pb)/V% have variances 

V(q^) = V(p + )= f r 1 . (95) 

The simplest description of the entanglement swapping 
protocol can be given when we consider the limit for /1 — > 
00. In this case the initial states are ideal EPR states 
with quadratures perfectly correlated, i.e., q a = qA and 
p a = —f>A for Alice, and q — qs and pb = ~pb for 
Bob. Then, the overall action of Charlie, i.e., the Bell 
measurement plus classical communication, corresponds 
to create a remote state with 

qb = q a ~ V2q-, p b = -p a - V%P+ ■ (96) 

The quadratures of the two remote modes are perfectly 
correlated, up to an erasable displacement. In other 
words, the ideal EPR correlations have been swapped 
from the initial states to the final conditional state p a b\-y- 



B. Entanglement swapping in a correlated 
Gaussian environment 

The theory of entanglement swapping can be extended 
to include the presence of loss and correlated noise. We 
consider our model of correlated Gaussian environment 
with transmission r, thermal noise ui and correlations G. 
The modified scenario is depicted in Fig. [TT] 
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for r > 1/2 and (1 — t) _1 < g < \g$\ there are separable 
AC environments which do the job. Thus, despite the 
presence of one-mode EB, the injection of separable cor- 
relations from the environment reactivates the swapping 
protocol generating remote EPR correlations. 

Finally, we have that SC environments (G —gl) never 
satisfy the EPR condition, so that the initial EPR corre- 
lations are always lost. However, as we shall see below, 
quantum entanglement can be swapped with success. 



FIG. 11: Entanglement swapping in the presence of loss, ther- 
mal noise and environmental correlations (correlated Gaus- 
sian environment). The Bell detector has been simplified. 



ng of EPR correlations 



For simplicity, we start by studying the evolution of 
the EPR correlations under ideal input conditions (fi — > 
+oo). The quadratures of the travelling modes, A and B, 
are transformed according to the the same input-output 
relations specified in Eqs. ((87)) and (|88|) . Then, the out- 
put modes A' and B' are measured by Charlie, while 
modes E[ and E' 2 are traced out. After the classical 
communication of the outcome 7, the quadratures of the 
remote modes a and b satisfy the asymptotic relations 



% = q a - \ - (q- - Vi - rS q , 



Pb = ~Pa -\I-[P+- Vi - tS p ) ) 



(97) 
(98) 



where S q = {qE t - qE 2 )/V% and 5 P = {pe^ +Pe 2 )/V2 are 
noise variables introduced by the environment. 

Considering the EPR quadratures q- and p+ , and av- 
eraging them over the states, we compute the asymptotic 
EPR variances 



V(q-) 



1 — T 



(col - ZG) 



Assuming the EB condition u) = webi we finally get 

A 00 , EB = i[(l+T)I-(l-r)ZG], 
T 



(99) 
(100) 



which is equal to Eq. (|92|) up to a factor t _1 . 

In the case of a memory less environment (G = 0) we 
get Aq^eb = (1 + T_1 )I > I; which means that the 
EPR correlations cannot be swapped to the remote sys- 
tems. However, it is evident from Eq. (|100p that there 
are choices for the correlation block G such that the EPR 
condition A^eb < I is satisfied. For instance, this hap- 
pens when we consider the AC environment (G =gZ). In 
this case it is easy to check that A^eb < I is satisfied 
for t > 1/4 and 



(1-r) 



< 9 < |5epr| 



(101) 



Under these conditions, the original EPR correlations are 
successfully swapped to the remote modes. In particular, 



2. Swapping and distillation of entanglement 

Here we discuss in detail how quantum entanglement 
can be distributed by the swapping protocol in the pres- 
ence of a correlated Gaussian environment. In particular, 
we aim to address the following questions: 

• Suppose that Alice and Bob are not able to share 
entanglement with Charlie because the environ- 
ment is one-mode EB, is it still possible for Charlie 
to distribute entanglement to the remote parties by 
exploiting the environmental correlations injected 
in the swapping protocol? 

• In particular, is the swapping successful when the 
environmental correlations are separable? 

• Finally, can Alice and Bob distill the swapped en- 
tanglement by one-way protocols? 

Our previous discussion on EPR correlations suggests 
that these questions have positive answers. Here we ex- 
plicitly show this is specifically true for quantum entan- 
glement by finding the corresponding good regimes of 
parameters for the Gaussian environment. 

In order to study the propagation of entanglement we 
first need to derive the CM V a &| 7 of the conditional re- 
mote state p a b\y. As before, we have two identical EPR 
states at Alice's and Bob's stations with CM V(/i) given 
in Eq. (|2"9"|) . The travelling modes A and B are sent to 
Charlie through a Gaussian environment with transmis- 
sivity r, thermal noise u> and correlations G. After the 
Bell measurement and the classical communication of the 
result 7, the conditional remote state at Alice's and Bob's 
stations is Gaussian with CM 



v a ^| 7 



l'l 



(M 2 - l)r 



( 



j_ 1 

1 e ' 1 
e 1 v 



1 A 

6' 



(102) 



where 



TH+(l-T){u-g), 



T/J,+ (l-T)(u+g') . (103) 



(See Appendix lAl for its derivation). 

From the CM of Eq. (TT02|) we derive the smallest PTS 
eigenvalue e quantifying the remote entanglement at Al- 
ice's and Bob's stations. In the limit of large input en- 
tanglement fi 3> 1, we find a closed formula in terms of 
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the environmental parameters, i.e., 
1 - T 



+ := e(r, u, g, g') , (104) 



which is equal to Eq. (f72j) up to a factor r^ 1 . As before, 
this eigenvalue not only determines the log-negativity but 
also the coherent information associated with the remote 
state Pab\~/- In fact, in the limit of large p, the determi- 
nant of the CM dTQU) becomes 



det V af ,| 7 = 



(l-r)/x 



{uj-g){Lo + g') = {e^ 2 . (105) 



Then, the reduced CM of Bob 

(,; 2 - l)r 



V, 



6| 7 



(106) 



has asymptotic determinant equal to det V^u = /i 2 . One 
can easily check that the symplectic spectra of V ob | 7 and 
Vf,| 7 are diverging for /i — > oo. As a result, we can use 
the formula of Eq. ([28} for the asymptotic coherent in- 
formation, which gives 



I(a)b) = log 



det V, 



b| 7 



e \ det V, 



loe 



a6|7 



ee 



(107) 



Thus, the PST eigenvalue of Eq. (|104l) contains all the 
information about the distribution and distillation of en- 
tanglement in the swapping scenario. For e < 1 entangle- 
ment is successfully distributed by the swapping protocol 
(log-negativity £ > 0). Then, for the stronger condi- 
tion e < 2e~ 1 ~ 0.73, the swapped entanglement can 
also be distilled into I(a)b) entanglement bits per copy 
by means of one-way protocols. Note that the condition 
for entanglement distillation for the indirect distribution 
(e < 2e~ 1 ) is weaker than that found for the direct dis- 
tribution (equal to e < e _1 ). 

Now, let us assume the condition of one-mode EB 
(lj = ljeb) so that the bipartite states before measure- 
ment p a A' and pB'b are separable (see Fig. [TTj) . We in- 
vestigate the entanglement in the remote modes a and b 
by computing the eigenvalue e(r, web, g, <?')■ I R the stan- 
dard memoryless case (G = 0) we have e = 1 + t _1 
which means that no entanglement can be swapped, as 
expected. For studying the general correlated environ- 
ment, we consider different numerical values of the trans- 
missivity r, and we plot the e(r, ueb, 3, <?') on the correla- 
tion plan. The results are shown in Fig. ll2l and are similar 
to those achieved for direct distribution (see Fig. [5]). 

In each panel of Fig. [T2l the bona-fide values for the 
correlation parameters are individuated by the non-black 
area. Remote entanglement is distributed (e < 1) for val- 
ues of the correlation parameters (g, g 1 ) belonging to the 
gray activation area. For r = 1/2, the activation area is 
populated by entangled environments only. The property 
that entangled environments are necessary for swapping 
entanglement is valid for any r < 1/2. In fact, suppose 
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FIG. 12: Analysis of the swapped entanglement e in the corre- 
lation plan for different values of the transmissivity r = 0.5, 
0.6, 0.75, and 0.9 (from top left to bottom right). Corre- 
sponding values of the thermal noise are determined by the 
one- mode EB condition uj = ojeb- In each inset, the non-black 
area identifies the set of physical environments, which are di- 
vided into separable (S) and entangled (E) environments. The 
gray region identifies the activation area, i.e., those environ- 
ments whose correlations are strong enough to generate re- 
mote entanglement (e < 1). Within the gray activation area, 
the environments that lie below the dashed curve are those 
able to generate distillable entanglement (e < 2e~ 1 ). 



that £ < 1 holds. By using its formula in Eq. (I104p and 
the bona-fide conditions of Eqs. (j4l)|) and (14*51) . we can 
write e 2 < 1 as 



- m' ' + u{g' - g) < 



(108) 



Now, for r < 1/2, we have r 2 (l — r) 2 < 1 and using 
this inequality in Eq. (|108[) . we derive 



gg' - 1< u{g - g') < u\g - g'\ 



(109) 



which is the entanglement condition for the environment 
[i.e., the violation of Eq. (|52"1) ]. 

It is clear that the most interesting result holds for 
transmissivities r > 1/2. In this regime, in fact, the 
distribution of remote entanglement can be activated by 
separable environments. As explicitly shown for r = 0.6, 
0.75 and 0.9, the activation area is progressively invaded 
by separable environments, with the entangled environ- 
ments almost disappearing for r = 0.9. In other words, 
separable correlations become more and more impor- 
tant for increasing transmissivities. Furthermore, for 
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r > 0.75, separable environments are even able to dis- 
tribute distillable entanglement (e < 2e _1 ). 

By comparing Fig. [5] and Fig. |T2l e.g., the two cor- 
responding insets for t = 0.75, it is evident how entan- 
glement is more easily generated by the direct protocol, 
as a consequence of the extra factor r" 1 in Eq. (|104[) . 
However, it is also evident that distillable entanglement 
is more easily generated by the swapping protocol, as 
consequence of the factor 2 in Eq. (|107j) which leads to 
the weaker condition e < 2e~ 1 . This is particularly ev- 
ident at high transmissivities (where r _1 becomes less 
important). 

It is important to note that our findings suggest a 
preferred working mechanism for quantum repeaters in 
the presence of correlated-noise environments, where the 
swapping protocol is performed before the entanglement 
distillation stage. In other words, Charlie first performs 
the swapping protocol on each pair of systems received 
from Alice and Bob. Correspondingly, the remote sys- 
tems at Alice's and Bob's stations are locally stored in 
quantum memories. Finally, these memories are manip- 
ulated by the coherent operations of the entanglement 
distillation protocol. This approach may be successful 
contrarily to the other strategy, based on distillation fol- 
lowed by swapping, which clearly cannot work in the 
present scenario. 

To conclude, it is easy to specialize the previous anal- 
ysis of the swapping protocol to the case of SC and AC 
environments, finding results similar to those found for 
the direct distribution. For brevity, here we consider only 
the extremal cases of the MSC and EPR environments 
(see Fig. [3]). As we know the MSC environments are 
optimal among the separable SC environments. Assum- 
ing one- mode EB, the remote entanglement generated by 
this environment is quantified by 



£msc 



V(l-r)(l + 3r) . 



(110) 



One can easily check that e < 1 if and only if 

r > (1 + \/5)/4 « 0.809 . (Ill) 

Thus, for sufficiently high transmissivities, the separable 
correlations of the MSC environment are able to activate 
the swapping of entanglement, despite the single-mode 
lossy channels are EB. It is also easy to check that dis- 
tillable entanglement (e < 2/e) is generated for 



e(e + 2VT + e2) n 
T> 4 + 3^ *°" 



(112) 



which is a regime easily accessible to Alice and Bob. 

For the positive EPR environment (g — |<7epr|) we 
have 



£epr. 



(113) 



which is less than 1 when r > 1/4. As expected, remote 
entanglement can be generated at transmissivities which 



are sensibly lower than those necessary for the MSC en- 
vironments. Furthermore, we can perform entanglement 
distillation (e < 2/e) for transmissivities 



r > 



(e-2)2 



0.289 



(114) 



VI. CONCLUSION AND DISCUSSION 

In conclusion, we have investigated the distribution 
of entanglement in the presence of correlated-noise envi- 
ronments, considering the framework of continuous vari- 
able systems, in particular, bosonic systems and Gaus- 
sian states. We have introduced and fully characterized a 
model of correlated Gaussian environment which gener- 
alizes the standard model of Gaussian environment with 
losses and thermal noise. 

Considering this correlated-noise environment, we have 
analyzed scenarios of direct distribution and indirect dis- 
tribution, i.e., entanglement swapping. In both cases, we 
have assumed the condition of one-mode EB, meaning 
that the transmission of a single bosonic mode cannot 
distribute entanglement between the parties, neither be- 
tween Charlie and Alice nor between Charlie and Bob. 
Despite this, we have shown that the distribution of en- 
tanglement is still possible when we consider the trans- 
mission of two bosonic modes, for instance, from Charlie 
to Alice and Bob in the scheme of direct distribution. 

The success of the two-mode transmission relies on 
the fact that environmental correlations are injected into 
the quantum systems and they can be strong enough to 
restore the entanglement broken by the thermal noise. 
While this process of reactivation is understandable in 
the presence of an entangled environment, where the orig- 
inal systems' entanglement is assisted or even replaced by 
the environmental entanglement, it is clearly paradoxical 
in the case of a separable environment, where no injection 
of entanglement may take place. 

Surprisingly, the injection of the weaker separable cor- 
relations is sufficient to restore the entanglement distri- 
bution, as we have shown for wide regimes of parame- 
ters. Furthermore, the generated entanglement can even 
be distillable by means of one-way protocols. The fact 
that separability can be exploited to recover from entan- 
glement breaking is clearly a paradoxical behavior which 
poses fundamental questions on the intimate relations 
between local and nonlocal correlations. In a few words, 
why do separable correlations have entangling power? 

In order to give more analytical examples, we have 
also defined two subclasses of correlated Gaussian en- 
vironments, corresponding to symmetric or asymmetric 
correlations among the quadratures. Using these specific 
environments, we have analyzed the evolution of the EPR 
correlations and studied the distribution of entanglement 
in terms of purely classical correlations and quantum dis- 
cord. We have found that entanglement can be restored 
by separable correlations which are mostly classical, i.e., 
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with a negligible amount ol quantum discord. 

Besides providing the first treatment of the problem 
of entanglement breaking in the presence of correlated 
noise, our work also extends the theory of entanglement 
swapping to this scenario, suggesting a preferred mecha- 
nism for quantum repeaters, where the swapping protocol 
must anticipate the distillation stage. 

In terms of potential impact, our analysis shows new 
perspectives for entanglement distribution and distilla- 
tion in the presence of memory channels and correlated- 
noise environments, such as those arising from the non- 
Mar kovian dynamics of open quantum systems (Hoi . l5l| . 
Note that memory channels and non-Markovian environ- 
ments are present in a wide series of practical scenar- 
ios. For instance, they naturally arise in the context of 
spin chains [201 ] and micromasers 21] . Other important 
examples can be found in condensed matter, in particu- 
lar when we consider the dynamics of quantum dots in 
photonic crystals [22| . In the bosonic setting, memory 
Gaussian channels come into play when electromagnetic 
modes propagate through dispersive media, such as linear 
optical systems or free-space. In this case, correlations 
and memory effects are naturally introduced by diffrac- 
tion [23T - I251 ] . Finally, other examples of bosonic mem- 
ory channels can be found in the radiation propagation 
through atmospheric turbulence (26l - [28j . 



symplectic matrix 



S = I®I®S(t)®S(t) 7 , 



(A4) 



where the identity matrices Iffil act on the remote modes, 
a and b, the beam splitter matrix S(r) [given in Eq. (jlll) ] 
acts on modes A and Ei, while S(r) acts on modes E 2 
and B. The second beam-splitter matrix is transposed in 
order to have positive reflection from the environment. 

The output state of modes abA' E[E 2 B' after the ac- 
tion of the interferometer is a Gaussian state with zero 
mean and CM equal to 



VabA'E'E'B' — S V ' a bAE 1E2B O 



After simple algebra, we get 



V a6 Wi W 2 

VabA'E^B' = I Wf V A , E[ W 3 



(A5) 



(A6) 



where the blocks along the diagonal correspond to the 
reduced CMs V Qb = p(I © I), 



V Vl , = I ' ! I I , V 



! - - I I ,-T ) ■ v - I zl tl ) ' (A7) 



with 



Appendix A: Computation of the covariance matrix 
of Eq. (fT02jl 

Consider the scenario depicted in Fig. [TT] We have a 
total of 6 input modes: Alice's modes a and A, Bob's 
modes b and B, and Eve's modes E\ and E%. The global 
input state is in a tensor product 



Pa A <8> PbB ® PE t E 2 



(Al) 



where p a A = PbB = p is an input EPR state with CM 
given in Eq. (|29[) and pe x e 2 is an environmental zero- 
mean Gaussian state with CM given in Eq. (pt2"j) . The 
global input state is therefore a zero-mean Gaussian state 
with CM 



V 



aAbBE 1 ^2 



V(/i) © V(m) © V Ei e 2 («, g, g') ■ (A2) 



It is helpful to permute the modes so to have the or- 
dering abAEiE 2 B, where the upper-case modes are those 
transformed by the beam splitters. After reordering, the 
input CM has the explicit form 
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abAExE-iB 
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fjfZ 


pi 










wi 
G 
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fj/Z 







(A3) 



V n'Z nl ) 



where is the 2x2 zero matrix. Now the global ac- 
tion of the two beam splitters can be represented by the 



x := rp + (1 - r)w , x' := tlj + (1- r)p , (A8) 



and z := y/r(l — t)(oj — p). The off-diagonal blocks are 
given by 



Wi 



W 2 = 



pl^Z -p'VT^ ~Z 











-p'y/T^ ~Z fi'yfFZ 



and 



W, 



Vr(l-r)G (l-r)G 

TG y/T{l-T)G 



(A9) 



(A10) 



(All) 



Since we are interested in the CM of Alice and Bob, 
we trace out the two environmental modes E[ and E' 2l 
which corresponds to deleting the corresponding rows 
and columns in the CM of Eq. (|A6I) . As a result, we 
get the following reduced CM for the modes abA'B' 



bA'B' 



pi p'y^FZ 

pi p'Vr~Z 

p'y/TZ xl (l-r)G 

p'y/r~Z (l-r)G xl 

( Al2 ) 

From this matrix, it is easy to check that the bipar- 
tite states before measurement p aJi > and pi,B> have CMs 
VaA' = Vt,B' of the same form given in Eq. It is 

also trivial to check that Alice and Bob's state pa'B' is 
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asymptotically separable for any r > 0. In fact, from the 
reduced CM ~Va'B'i we can compute the smallest PTS 
eigenvalue 

e A >B> =T l x+^^{2u-\g-g'\) , (A13) 

which is always > 1 for large /i and r > 0. In the singular 
case t = 0, we have pa'B' — Pe ± e 2 so that entanglement 
or separability directly comes from the environment. 

Most importantly, from the CM of Eq. (|AF2l) we can 
compute the CM V ab | 7 of the conditional remote state 
p a b| 7 after the Bell measurement. In order to achieve this 
result, we use the transformation formula for CMs under 
Bell measurements which has been proven in Ref. [52j . 
As a first step, we put "V a bA' B' in the blockform 

f A d C 2 \ 
VabA>B> =\CJB 1 D , (A14) 
V D r B 2 ) 

where 

A = (jl pi, Bi = B 2 = xl , D = (1 - r)G , (A15) 
and 




(A16) 



From the blocks of the CM (|A14I) . we construct the fol- 
lowing matrix 

= i(ZBiZ + B 2 ZD D T Z) , (A17) 

where Z is the reflection matrix of Eq. (130)) . Then the 
conditional CM is given by the formula (52| 

Vab ^ = A _ 2det0 £ Ci(Xf 0X,-)Cj , (A18) 
where 

Xi = ( j 1 ) , X 2 := ( _ x 1 ) = O . (A19) 
Explicitly, we compute the theta matrix 

(A20) 

with elements given in Eq. (|103[) . Using this matrix in 
Eq. (|A18|) we get the final expression of Eq. (|102[) . It is 
easy to check that for r = 1 (absence of loss and noise) , 
we get Eq. (1941) which is the result already known in the 
literature [49i. l53l|. 
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